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Cn ' Abstract - In this paper we investigate a tachyon field model in cosmology, provided its 

interaction with the quintessence or phantom fields. The model takes into account this inter- 
action beyond the usual approach, in which the interaction is phenomenologically described 
by the energy flow between the matter components. In our model, the interaction of tachyon 
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jH ■ field with a canonical scalar field is taken into account through the interaction potential in 

^^ ' the total Lagrangian of the system, like in the case of two or more canonical scalar fields. We 

obtain the different types of exact solution for the model by employing the so-called "first 
_-•. ■ order formalism" procedures. 

(N 

PACS: 98.80.-k,98.80.Jk. 
O ' Key words: Cosmological model, tachyon field, interaction, exact solutions. 

u 



> 



(N 



X 



1 Introduction 



Several cosmological observations, such as SNe type la [J [5], Cosmic Microwave Background 

f^ , Radiation [3J H], Baryon Acoustic Oscillations in galaxy surveys [5J H] etc., indicate that the 

(T^ ' observable Universe is undergoing a phase of accelerated expansion. In order to explain such 

^D . phenomena, one needs to add an extra energy component of the Universe, which has a large negative 

^'J ' pressure and also dominates over all other components e.g. non-relativistic matter, radiation etc. 

C^ . As known, this mysterious component is usually referred to as "dark energy" (DE). 

C^ ■ To accelerate the Universe expansion, equation of state (EoS) of DE, w = p/p where p and p 

are the pressure and energy density respectively, must satisfy w < —1/3. The simplest candidate 

for the role of DE is the cosmological constant, for which EoS is w = — 1. However, there is some 

evidence that DE may evolve from w > —1 in the remote past to w < — 1 at present. To date, there 

have been studied a large class of scalar - field models of DE [7] , including tachyon [5] , phantom 

^ , [5], quintom [TU] and other models. Furthermore, the proposals for DE include interacting DE 

^_' model [TT], braneworld models [H] and holographic DE mo dels P^. Phantom field (with negative 

kinetic energy) was also supposed as a candidate for DE, since it allows sufficient negative pressure 

with EoS w < —1 [m [TS]. The so-called quintom scenario of DE has been developed in order to 

explain the dark energy EoS crossing —1. The first model of quintom scenario of dark energy has 

been obtained in [TU]. This model is studied in detail later e.g in [TBI [TT] . 

Recently there has been increasing interest for constructing the tachyon models in cosmology 
[18] , where the appearance of a tachyon field is mainly motivated by the string theory [19] . In the 
context of string theory, the tachyon field in the world volume theory of the open string stretched 
between a D-brane and an anti-D-brane or a non-BPS D-brane plays the role of scalar field in 
the quintom model. The effective action used in the study of tachyon cosmology consists of the 
standard Einstein-Hilbert action and an effective action for the tachyon field on unstable D-brane 
or D-brane anti D-brane system. What distinguishes the tachyon action from the standard Klein- 
Gordon form for scalar field is that the tachyon action is of the " Dirac-Born-Infeld " (DBI) form. 
As known, tachyon fields are able to contribute for the inflation and late time acceleration |20) . 
The tachyon model of inflation has been discussed in [5T1 [22j . Tachyon dark energy as a function of 
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various forms of its potential was studied in [13]"I2S]- A model for tachyon field with non-minimal 
derivative coupling to gravity has been considered in [27] . 

In view of the latest developments in DE modeling, the study of interacting tachyon models 
acquired much significance. It seems natural to study different hybrid models which could include 
an interacting tachyon field in the framework of different approaches. However, the model ac- 
counting for the interaction of tachyon field with other fields has a well-known problem. So far the 
description of tachyonic interaction is made phenomenologically by involving an energy exchange 
term to govern the interaction (see, e.g., p8 ] -[44 ] ). 

In this paper we made an attempt to develop an alternative approach for description of the 
interacting tachyon field, which is based on a more conventional phenomenological description 
via an interaction potential. We obtain the different types of exact solution for the Friedmann- 
Robertson- Walker (FRW) model filled with the interacting tachyon and quintessence (phantom) 
scalar fields by employing the so-called first order formalism procedure |451 146] . 



2 Basic equations 

An effective DBI Lagrangian density inspired from the super-string theory is 

CDBi = -r\^)4^^~Jmy^+r\4>)-ui{4>), (i) 

where the tension of the brane is /^^(0) — T^h'^{(t>) with /i(0) being the warped factor in the the 
metric 

dsla = h'^{r)dsl + h-^{r){dr^ + r^ds^), (2) 

and Ui{(j)) is the potential for the DBI-field, arising from quantum interactions between the D3- 
brane associated with 0, and other D-branes. In terms of the brane tension T((/>), the Lagrangian 
density is 



L^ = ~r(0)^I - (f^/Ticj,) + T{4>) - Uii^), (3) 

for a homogeneous and isotropic Universe. By adding to the tachyon component a scalar field tp 
described by the Lagrangian density with an interaction potential 

we have the total Lagrangian density C = C^ + C^ for the hybrid source of gravity. Here, e — +\ 
represents quintessence while e — —\ refers to phantom field. 

Relying on a phenomenological approach in the presence of scalar field, we assume further that 
the brane tension depends also on (/? as a parameter. If we substitute r(0, ip) — Ui{(j>) + U2{p) + 
Uint{4>i v) iiito the total Lagrangian density and denote T{(f), ip) — V{(J3, ip), i.e. 

y(0,^) - t/i(0) + U2{p) + U„,t{^,p), (5) 

then we can write down the following Lagrangian density for our model: 



Some motivation for this model is that in the case of slow roll, viz. ip^ << V{(l),p), the total 
Lagrangian is reduced to the canonical quintom model. 

Keeping in mind that all quantities, 4>{t),p{t) and V {(j){t) ^ p{t)) are functions of time t which 
should be found from the set of field equations, we can suppose the re-definition of tachyon field 
as follows (cf [47): 

Substitution of the latter into © leads to the following Lagrangian density: 

C = -V{x,p)Vl^+'^^'- (8) 



For the sake of simplicity, we consider a spatially flat FRW cosmological model with the space-time 
interval 

ds^ = dt^ - a^{t) (dr^ + r'^dO?) , (9) 

where a{t) is a scale factor. Taking into account this metric, we can derive from Lagrangian ([5]) 
the following equations for the tachyon and scalar fields: 



X 



i~x' 



3iJx+^+XV^^=0, 



(^ + 3iI^ + eF^A/l-X =0, 



(10) 



(11) 



where H = d{t)/a{t) is the Hubble function, V-^ — dV[xi '^)/dXi and V^ = dV[xi '^)/dv- Besides, 
with the help of well-known formula for the energy-momentum tensor, T^^ = 2dC/dg^'^ — 9fj.u^, 
and from ([8| it follows that the total effective energy density and pressure are 



ptot = ^^= + ^^', Vtot = -v{x, v)^i-x' + -^' 



(12) 



Hence, the behavior of this model is defined by equations (fTOl) . (fTTj) and the Friedmann equations, 
(3/2)iJ^ = Ptot and H ~ —{ptot + Ptot) where AttG = 1. From the metric ([9]) and equations (|12p . 
one can conclude that 



H' 



v{x,v) 



x' 



-.V 



v/i-x' 

At the same time, the energy-momentum conservation law. 



(13) 
(14) 



Ptot + 3H{ptot + Ptot) = 0, 



leads to the following equation 



xv 



VT 



X 



1-r 



+ 3i?x- 



V 



•^ 



e{(p + iHiP) 



K. 



V^^' 



= 0, 



(15) 



which is satisfied due to the field equations (ITUl) . ((TT|) . 



3 Exact solutions 



3.1 The method 

Even for the given potential V{x, v), it is difficult to find exact solutions for our model. However, 
a class of exact solutions can be obtained in terms of the so-called superpotential in the first order 
formalism. This procedure was first performed with a single scalar field in |45) . and it was later 
re-opened and extended on the case of two or more fields with standard dynamics in [46]. The 
method of superpotential can be also productively applied to the quintom models (see, e.g., [IS] 
and [in]). 

It should be noted that the general first-order equations for multiple fields, that are covariant 
in field space, were written down in [50| . As seems, there exist some examples for which the first- 
order framework, as written in terms of a superpotential, fails [51]. Nevertheless, we are able to 
give several examples for which the method works. 

Here, we introduce the superpotential function W{x, f) by the equality 



H = W{x,^), 



(16) 



in which the Hubble parameter H(t), as a function of time, is presumably expressed in terms of 
fields x{'t),'p{t)- By inserting (IT51) into (|14l) . one can obtain two first-order equations as follows: 



2W^ 



3W^ 



eW^ 



, ip = -eW, 



v 



(17) 



where W^ = dW/dx, W^ = dW/dip. The potential y(x, (^) is followed from HH), (dH) and (H?]) 
in the form: 

V ^ y {3W' - eW^y - iW^ (18) 

The evolutionary equation for the superpotential W(x{t),(p{t)) is followed from W = W^x + 
W^ip and Eq. (HH) as 

VF = ^ eW"^. (19) 

With the help of Eqs. (O, P7|) and (IT51) . it is easy to find that the original tachyon field (/)(i) can 
be obtained from the following equation 



AW?, 



?=2W^Jl- ,,„,, \,,,,, = 2W^ Vl-X"- (20) 



(31^2 _ g^2)2 



Ptot _ ^ 2 H 

Plot 3 ^2 

H 



To analyze this model, it is useful to consider the EoS parameter w — = — 1 ^ and 



the deceleration parameter q = —^^ = —1 — — r- followed from Eqs. (fldl) . ([19)) : 



{:}--{r} 



W^2f3Pl/2_£W2^ 14^2 



(21) 



So we have a wide range of possibilities to solve the model equations assuming some certain 
dependence W{xt'p)- Instead, we can provide several classes of solution for the model evolving 
from some conditions on superpotential. Below, we show how it can be realized with the help of 
some ansatz for the superpotential. 

3.2 The first ansatz 

The simplest condition on superpotential may be represented as 

W^^X^W'^ix,^), W^ = \2W'^{x,^), (22) 

where m, Ai, A2 are constants. Then from (fTTll and ([22]) . we have 

0\ ^iJim — 2 

In view of Eqs. (|19p and (|22p . we can obtain the evolutionary equation for the superpotential or, 
keeping in mind Eq. ([TS]). the Hubble parameter as follows 

2 T;f/2m-2 



2\iW 



W = - 3_;a2vK2"^-2 - ^^2 W^-. H = W{t). (24) 



Integrating this equation for a given to, we can find the fields from Eq. (1^^ . Then it is possible to 
find the field potential (|18p . the Eos index and the deceleration parameter (pij) . Let us illustrate 
this procedure for two different values of to. 

3.2.1 The case m = 

In this case, it follows from Eq. ([22]) that W^ ~ Xi,W^ ~ A2. Therefore, 

Wix,ip) = XiX + ^2^, (25) 

and the field equations (|23|) become as follows 



From Eq. ([M]) . one can find the following equation for the Hubble parameter (or superpotential) 

2X1 



H 



3i?2 _ £A2 



eA^ 



which can be integrated in an implicit form as 

1 „ 2\\ 



t_ ' H \ -■' / *""''_ 1 ( 3A2g ■ 

Ai A3^3A|-6A^l tanh ^ J V^/3Af— 6Af, 



for e = +1 and e = — 1, consequently. Besides, we can obtain that 



{:}--{r} 



2A2 



>2 
^^2 



H'^iZm-eXl) m 



(27) 



(28) 



(29) 



3.2.2 The case m = 1 

From Eq. ([22]) with m = 1 and dW = VK^^X + VK;pC?i^ , it follows that 

W{x, f) = Wo exp(Aix + A2V3), 
Then from ([23l) . we have 



X = 



2Ai 
'3~^ 



W^\ ip = -eX2W. 



With the help of Eqs. ([T8| and ([22]) . we can find the following expression for the potential: 



W 



V=—^{i-eXiYW'-AXl 



(30) 
(31) 

(32) 



where W is presented by (|30|) . Besides, we are able to obtain the equation for the superpotential 
followed from ([21]) : 

2A2 



W = -{tXi,W^ 



'i-eXl 



Integrating Eq. p3p with e = — 1, we have 

H ^ Ho coth \xlHo{tf - t) , 



(33) 



(34) 



A / 2 

where Hq — -r-\ rj? ^-nd t/ > is a constant of integration. Therefore, the scale factor can 

A2 V 3 + A, 



be written down as 

a{t) = aol sinh 
From Eqs. dH]) and dMl), we can find that 



A^i/o(t/ - t)] } 



Ai 



cosh 



XrjHl 



o{tf-t)]} 



(35) 



(36) 



In view of Eqs. (|3ip and (|34| . it can be found that 



Aix = In cosh X\Ho{tf — t) 



X2(p — In sinh 



\Ho{tj — t) 



+ c, 



(37) 



where C is a constant. The behavior of this model in time is plotted in Fig. 1. 



H 




Figure 1: The plot of the scale factor a, the Hubble parameter H, the EoS parameter w, and the 
deceleration parameter q versus time t for to = 1, e = —1, Ai = \f2, A2 = 1, i/ = 2, and ao = 2. 



3.3 The second ansatz 

Let us now consider the case when 

In view of this ansatz and equation (J17p , we have 

2X' $ 



X- 



X^ 3$2-e$'2' 



■^ 



-eX$', 



where the primes denote derivatives with respect to the appropriate arguments. 
In order to obtain an exact solution, we suppose that 

W = \^, 

X 

where A is a constant. Hence, we have the following equation for the potential: 

^(X: ^) = ^ VA^ (3(^2 _,)2_ 4^2. 

From Eqs. (l39l) and (|40)) . we can find the following equations: 



X 



2lp 



A 



A(3^2 „ g) 
This set of equations can be written down as follows 



<^ 



(p = B 3^2-e 



l/SeA^ 



X 



V 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



where B ^ is a constant of integration. The first equation in (j45|) is exactly differentiable for 
several particular values of A. Let us consider two simple examples of exact solution for e = ±1. 

3.3.1 The case e = +1, A = l/%/3 

Integrating the first equation in (j43p with x > 0, we can obtain that 



Lp = —j= tanh VSBitf ~ t) 
V3 L 



, X 



cosh 



V3B 



V3B{tf - t) 



(44) 




Figure 2: The plot of the scale factor a, the Hubble parameter H, the EoS parameter w, and the 
deceleration parameter q versus HqI for e = — 1, A = l/-\/3, and Hq = 18. 



Due to Eqs. (gO]) and (gH), we have 



Hit) 



B s\nh[^/iB{tf - t)] 



a{t) = flo exp ■ 



VS cosh!^ [V3B{tf -t)]' 
1 



6cosh^[V5B{tf -t)] j ' 
Finally, we can obtain the following expressions for the EoS and deceleration parameters: 

{™} = -l-{g}x cosh^[V3B{tf - i)]| coth^[V3B{tf - t)] - s}. 



(45) 
(46) 

(47) 



3.3.2 The case e= -1,A = 1/V3 

In this case, we can integrate equation (j43p to obtain 

1/3 



Viip = (Hat + -^m^) - (Hot + ^WHlt'^^ 



-1/3 



(48) 



where, for the sake of simplicity, we let the integration constant be zero, and put B = 2Hq/3\/3 
with Ho being a constant. Then it is easy to find that x = 1/ v3<P and 



Hit) 



Ho 



{ (Hot + ^l + H^) ^"' - (Hot + ^1 + H^ot") '^'} 



9^1 + H^t^ 
By using (|2T|) and p9)) . we are able to find the EoS and deceleration parameters as follows: 



p. _ r2/3| f^^^ 



6g^+Ht) 
2g(-)it) [gi-Ht)]^ 



where 



or 



2/3 



ff(±)(t) = (Hot + ^1 + H§t^y'"± (Hot + ^l + H^f^ 



-2/3 



(49) 



(50) 



(51) 



By using the obvious equality 

it is easy to integrate Eq. (|49l) and obtain the following result for the scale factor: 

a(i)=aoexp[ig(+)(Ol, (52) 

L6 J 

where oq is a constant of integration. The main features of this solution is shown in Fig. 2. It is 
interesting to note that, passing through a maximum, then iJ — > 0, and w, g — > —1 as i ^^ oo. 

4 Conclusion 

In this paper, we have studied the flat FRW cosmological models with interacting quintessence 
(phantom) and tachyon scalar fields considered as the origin of gravity. We have described theo- 
retical models based on the assumption that the interaction of tachyon field could be described with 
the help of potential. We have obtained the different types of exact solution for the Friedmann- 
Robertson- Walker model filled with the interacting tachyon and quintessence (phantom) scalar 
fields by employing the so-called first order formalism procedure. 

We have left aside the question of the possibility to reconstruct the potential in the form of an 
explicit function of the original tachyon field (j){t). This issue will be studied further. 

Of course, all given examples are not the limit of the capacities of the method considered. We 
hope that the derived model is the next step in the development of tachyon cosmology, and it can 
be utihzed to describe the dynamics of the evolution of the actual Universe. 
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